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REDUCED COMPLEXITY CHANNEL ESTIMATION FOR WIRELESS 

COMMUNICATION SYSTEMS 

CROSS-REFERENCE TO RELATED APPLICATION 
[0001] This application claims the benefit of U.S. Provisional Patent Application Serial No. 

60/427,896, filed November 19, 2002, which are incorporated herein by reference in its entirety. 

BACKGROUND 

/. Field 

[0002] The present invention relates generally to data conmiunication, and more specifically to 

techniques for performing channel estimation with reduced complexity. 

//. Background 

[0003] Wireless communication systems are widely deployed to provide various types of 

communication such as voice, packet data, and so on. These systems may be multiple-access 
systems capable of supporting communication with multiple users by sharing the available 
system resources. Examples of such multiple-access systems include code division multiple 
access (CDMA) systems, time division multiple access (TDMA) systenis, and orthogonal 
frequency division multiple access (OFDMA) systems. 

[0004] Orthogonal frequency division multiplexing (OFDM) effectively partitions the overall 

system bandwidth into a number of (N) orthogonal subbands. These subbands are also referred 
to as tones, frequency bins, and frequency subchannels. With OFDM, each subband is associated 
with a respective subcarrier upon which data may be modulated- Each subband may thus be 
viewed as an independent transmission channel that may be used to transmit data. 

. — [0005]- In-a wireless -communication system,-an RF modulated signal-from a transmitter may 

reach a receiver via a number of propagation paths. For an OFDM system, the N subbands may 
experience different effective channels due to different effects of fading and multipath and may 
consequently be associated with different complex channel gains. 

[0006] An accurate estimate of the response of the wireless channel between the transmitter and 

the receiver is normally needed in order to effectively transmit data on the available subbands. 
Channel estimation is typically performed by sending a pilot from the transmitter and measuring 
the pilot at the receiver. Since the pilot is made up of symbols that are known a priori by the 
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receiver, the channel response can be estimated as the ratio of the received pilot symbol over the 
transmitted pilot symbol for each subband used for pilot transmission. 

[0007] Pilot transmission represents overhead in a wireless communication system. Thus, it is 

desirable to minimize pilot transmission to the extent possible. However, because of noise and 
other artifacts in the wireless channel, a sufficient amount of pilot needs to be transmitted in 
order for the receiver to obtain a reasonably accurate estimate of the channel response. 
Moreover, the pilot transmission needs to be repeated to account for variations in the channel 
over time due to fading and changes in the multipath constituents. Consequently, channel 
estimation normally consumes a noticeable portion of the system resources. 

[0008] In an OHDM system, to reduce the amount of overhead for pilot, a pilot transmission may 

be sent on a group of designated subbands, which may be only a subset of the available 
subbands. An initial estimate of the channel response may be obtained for the designated 
subbands based on the pilot transmission. Signal processing may then be performed to obtain an 
enhanced channel response for a group of desired subbands, which typically includes the 
subbands to be used for data transmission. The signal processing may further perform noise 
averaging to obtain a more accurate estimate of the channel response. As described in detail 
below, depending on the number of designated subbands used for pilot transmission and the 
impulse response of the channel, the signal processing may be computationally intensive and 
require a large number of complex multiplications. 

[0009] There is therefore a need in the art for techniques to more efficiently derive an estimate of 

the channel response in a wireless communication system, such as an OFDM system. 

SUMMARY 

[0010] Techniques are provided herein to derive a channel estimate using substantially fewer 

numbers of complex multiplications than with a brute-force method to derive the same channel 
estimater This channel estiniate niay be a leasT^^^ of the impulse response of a 

wireless channel, which may be derived based on an initial frequency response estimate H of the 
wireless channel. As described in detail below, the least square channel impulse response 
estimate may be derived by a matrix multiplication between the vector H and a matrix w'* , 
which is derived based on a discrete Fourier transform (DFT) matrix W . The structure of the 
matrix W can be exploited to decompose the matrix multiplication W H into a sum of matrix 
multiplications between smaller sub-matrices of W and smaller sub-vectors of H. The 
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properties of the sub-matrices of W can be exploited to simplify the computation. The net result 
is fewer number of complex multiplications required to obtain the least square channel impulse 
response estimate. 

[0011] In one embodiment, a method is provided for deriving an estimate of a wireless channel 

in a wireless communication system (e.g., an OFDM system). In accordance with the method, an 
intermediate vector B is initially obtained, which is derived based on K sub- vectors of the vector 

H for a first channel estimate (e.g., a channel frequency response estimate) and at least two DFT 
sub-matrices for the DFT matrix W , where K is an integer greater than one. An intermediate 
matrix A for the DFT matrix W is also obtained. A second channel estimate (e.g., a least 
square channel impulse response estimate) is then derived based on the intermediate vector B 
and the intermediate matrix A . 

[0012] In one implementation, the intermediate vector B is obtained by first computing DFTs of 

a first matrix H^j^l* which is formed based on the vector H, to provide a second matrix Gl^^l- 
Inner products between the columns of a base DFT sub-matrix and the rows of the second 
matrix Gl^l then computed to obtain the entries of the intermediate vector B . Details of this 
implementation are described below. 

[0013] Various aspects and embodiments of the invention are described in further detail below. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0014] The features, nature, and advantages of the present invention will become more apparent 

from the detailed description set forth below when taken in conjunction with the drawings in 
which like reference characters identify correspondingly throughout and wherein: 
[0015] FIG. 1 shows an OFDM subband structure; 

- - [0016] FIG. 2A shows the relatidhship between the response and the impulse response 

of a wireless channel; 
[0017] FIG. 2B graphically shows a DFT matrix W ; 

[0018] FIG. 3 A graphically shows the relationship between DFT matrices W and W ; 

[0019] FIG. 3B graphically shows the partitioning of the vector H into K sub-vectors and the 

DFT matrix W into K sub-matrices; 
[0020] FIG. 3C graphically shows the relationship between the matrices Wlx^ and H^xl • 
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[0021] FIGS. 4 and 5 show two processes for deriving the least square channel impulse response 

estimate using a low complexity least square method; and 
[0022] FIG. 6 shows a block diagram of an access point and a terminal in a wireless 

communication system. 

DETAILED DESCRIPTION 
[0023] The channel estimation techniques described herein may be used for any communication 

system with multiple subbands. For clarity, these techniques are described specifically for an 
OFDM system. 

[0024] FIG. 1 shows a subband structure 100 that may be used for an OFDM system. The 

OFDM system has an overall system bandwidth of W MHz, which is partitioned into N 
orthogonal subbands using OFDM. Each subband has a bandwidth of W/N MHz. In a typical 
OFDM system, only M of the N total subbands are used for data transmission, where M < N . 
These M usable subbands are also referred to as data subbands. The remaining N - M subbands 
are not used for data transmission and serve as guard subbands to allow the OFDM system to 
meet spectral mask requirements. The M usable subbands include subbands F through 
F + M-1. 

[0025] For OFDM, the data to be transmitted on each subband is first modulated (i.e., symbol 

mapped) using a particular modulation scheme selected for use for that subband. Signal values 
of zero are provided for the N-M unused subbands. For each symbol period, the M 
modulation symbols and N-M zeros for all N subbands are transformed to the time domain 
using an inverse fast Fourier transform (IFFT) to obtain a "transformed" symbol that includes N 
time-domain samples. The duration of each transformed symbol is inversely related to the 
bandwidth of each subband. For example, if the system bandwidth is W = 20 MHz and 
N = 256 , then the bandwidth of each subband is 78.125 KHz (or W/N MHz) and the duration 
of each transformed symbol is 12.8 jiisec (or N/ W fisec). 

[0026] The N subbands of the OFDM system may experience different channel conditions (e.g., 

different effects due to fading and multipath) and may be associated with different complex 
channel gains. An accurate estimate of the channel response is normally needed in order to 
properly process (e.g., decode and demodulate) data at the receiver. 

[0027] For clarity, in the following description, lower-case letters are used for indices, upper- 

case letters are used for constants, and bolded and underlined lower and upper-case letters are 
used for vectors and matrices. 
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[0028] The wireless channel in the OFDM system may be characterized by either a time-domain 

channel impulse response^ h , or a corresponding frequency-domain channel frequency response, 

H. The channel frequency response H is the discrete Fourier transform (DFT) of the channel 

impulse response h . This relationship may be expressed in matrix form, as follows: 

H = Wh , Eq(l) 
where h is an (Nxl) vector for the impulse response of the wireless channel between a 

transmitter and a receiver in the OFDM system; 
H is an (Nxl) vector for the frequency response of the wireless channel; and 
W is an (NxN) DFT matrix used to perform the DFT on the vector h to obtain the 

vector H. 

The DFT matrix W is defined such that the (n,m) -th entry, , is given as: 

(n-lKm-l) 

w^^=e N ,forn = {l ... N} and m = {l ... N}, Eq (2) 

where n is a row index and m is a colunm index. 

[0029] The vector h includes one non-zero entry for each tap of the channel impulse response. 

Thus, if the channel impulse response includes L taps, where L < N , then the first L entries of 
the vector h would be L non-zero values and the (N - L) following entries would be zeros. 

[0030] FIG. 2A graphically shows the relationship between the channel frequency response H 

and the channel impulse response h. The vector h includes N time-domain values for the 
impulse response of the wireless channel from the transmitter to the receiver, where some of the 
entries in h may be zeros. This vector h can be transformed to the frequency domain by pre- 
multiplying it with the matrix W . The vector H includes N frequency-domain values for the 
complex channel gains of the N subbands. 

[0031] FIG. 2B graphically shows the matrix W , which is an (N x N) matrix comprised of the 

.... elements for n = {1 ...-N}— and~m = {1 .7. N}-, which are defined in equation' (2). The 

superscript " " is not shown in FIG. 2B for clarity. Each row of the matrix W corresponds to 
one of the N total subbands. 
[0032] The impulse response of the wireless channel can be characterized by L taps, where L is 

typically much less than the number of total subbands (i.e., L< N). That is, if an impulse is 
applied to the wireless channel by the transmitter, then L time-domain samples (at the sample 
rate of W) would be sufficient to characterize the response of the wireless channel based on this 
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impulse stimulus. The number of taps L for the channel impulse response is dependent on the 
delay spread of the system, with a longer delay spread corresponding to a larger value for L. 

[0033] Because only L taps are needed for the channel impulse response, the channel frequency 

response H lies in a subspace of dimension L (instead of N). More specifically, the frequency 
response of the wireless channel may be fully characterized based on the channel gains for as few 
as L appropriately selected subbands, instead of all N subbands. Even if more than L channel 
gains are available, an enhanced estimate of the frequency response of the wireless channel may 
be obtained by suppressing the noise components outside this subspace. 

[0034] In one channel estimation technique, a more accurate estimate of the frequency response 

of a wireless channel is obtained based on a 3-step process. In the first step, an initial estimate of 
the channel frequency response,. H, is obtained based on the received and transmitted pilot 
symbols for each of S designated subbands, where S is an integer selected such that L < S < M . 
The S designated subbands may include all or only a subset of the M usable subbands. The 

initial channel frequency response estimate, H, may be expressed as: 

H = r,/x,=H,+n,/x, , Eq (3) 

where is a "receive" vector with S entries for the symbols received on the S designated 
subbands; 

X, is a "transmit" vector with S entries for the symbols transmitted on the S designated 
subbands; 

is an (Sxl) vectors that includes only S entries of the (Nxl) vector H for the S 
designated subbands; 

is a vector with S entries for additive white Gaussian noise (AWGN) received on the S 
designated subbands; and 

_ §5 /b^_= [q, /^i ^2 7^2 . • J^s /^sl^-'^^^^ includes -S ratios- for. the S- designated 

subbands. 

[0035] In the second step, a least square estimate of the impulse response of the wireless channel, 

h , is obtained based on the following optimization: 

h'=min||H-Wh,.||^ , Eq (4) 

where hy is an (Lxl) vector for a hypothesized impulse response of the channel; 
W is an (SxL) sub-matrix of the (NxN) matrix W; and 
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h is an (Lxl) vector for the least square channel impulse response estimate. 

[0036] FIG. 3A graphically shows the relationship between the matrices W and W . The S 

rows of the matrix W are the S rows of the matrix W corresponding to the S designated 

subbands. The L colunms of the matrix W are the first L columns of the matrix W . 
[0037] The solution to equation (4) that results in the minimum mean square error (or more 

specifically, the minimum Euclidean norm) may be expressed as: 

h'' = (W" W)-^ W" H = W*" H , Eq (5) 

where is an (LxS) matrix defined as w'' = (w'' W)"^ w" . 
[0038] In the third step, an enhanced estimate of the frequency response of the wireless channel, 

^ Is ^Is 

H , is obtained based on the least square channel impulse response estimate, h , as follows: 

^ Is ^Is 

H =Wh , Eq(6) 
where W is a (QxL) sub-matrix of the (NxN) matrix W; and 

h'^ is a (Qxl) vector for the enhanced channel frequency response estimate for Q 
desired subbands. 

The Q rows of the matrix W are the Q rows of the matrix W corresponding to the Q subbands 
for which the enhanced channel frequency response estimate is desired. In general, the matrix 
W may include any number and any combination of rows of the matrix W . For example, the 

matrix W may include only the S rows of the matrix W , the S rows of the matrix W plus one 
or more additional rows, the M rows of the matrix W for the M usable subbands, and so on. The 
group of S designated subbands may thus be the same or different from the group of Q desired 
^ subbands — - - 

A Is 

[0039] Equation (6) indicates that the enhanced channel frequency response estimate H may be 

obtained for Q desired subbands based on the least square channel impulse response estimate h 

that includes only L entries, where L is typically less than S and Q and may be much less than S 
and Q. 

[0040] The 3-step channel estimation technique is described in further detail in U.S. Patent 

Application Serial No. [Attorney Docket No. PD020718], entitled "Channel Estimation for 
OFDM Conmiunication Systems," filed October 29, 2002. The reduced complexity channel 
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estimation techniques described herein may also be used in conjunction with pilot transmission 
schemes described in U.S. Patent Application Serial No. [Attorney Docket No. PD020656], 
entitled "Uplink Pilot and Signaling Transmission in Wireless Conmiunication Systems," filed 
October 29, 2002. Both of these patent applications are assigned to the assignee of the present 
application and incorporated herein by reference. 
[0041] An OFDM system may be designed with a relatively large number of subbands. For 

example, an OFDM system may be designed with 256 total subbands (i.e., N = 256) and 224 
usable subbands (i.e., M = 224). In an example design, S may be selected to be equal to 224 for 
the downlink (i.e., S^, = 224 ) and equal to 32 for the uplink (i.e., = 32 ). The number of total 

subbands may be given with respect to L such that N = L*T. The number of designated 
subbands may also be given with respect to L such that S = K • L . For the example design 
described above with L = 16 and T = 16, K would be equal to 14 for the downlink and to 2 for 
the uplink (i.e., S^, = 14L and S„, = 2L). 

[0042] A straightforward or brute-force method for deriving the estimate h using equation (5) 

would require C^f =L-S complex multiplications for the matrix multiply between the (LxS) 

^ Is ^ ■ ^Is 

matrix W and the (Sxl) vector H. This is because each of the L elements of the vector h 
requires S complex multiplications for the inner product between one row of the matrix W and 
the vector H. For the example OFDM system described above, the number of complex 
multiplications required to derive - the estimate h can be given as 
= L • L • K = 16 • 16 • K = 256K , where K = 14 for the downlink and K = 2 for the uplink. A 

large number of complex multiplications may thus be required to derive the estimate h , 

especially for the downlink. _ _ 

[0043] Techniques are provided herein to derive the estimate h using substantially fewer 

numbers of complex multiplications than with the brute-force method. The structure of the 
matrix W can be exploited to decompose the matrix multiplication W H in equation (5) into a 
sum of K matrix multiplications between smaller sub-matrices of W and smaller sub-vectors of 
H. The properties of the sub-matrices of W can be exploited to simplify the computation. The 
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net result is fewer number of complex multiplications required to obtain the estimate h , as 
described in detail below. 

[0044] The (Sxl) vector H, where S = K-L, can be partitioned into K smaller (Lxl) sub- 

vectors, as follows: 



H = 



H, 



H. 



Eq(7) 



Each sub- vector , for A: = {1 ... K} , may be expressed as: 



Hp 



Eq(8) 



'F+(*-l)L ^^F+(*-l)L+l 

where H^^^f^_^^^^j is the estimated channel gain for subband F+(^-l)L+7, which may be 

obtained as shown in equation (3); 
F is the index for the first usable subband, as shown in FIG. 1; and 
" ^ " denotes the transpose. 

[0045] The (SxL) matrix W, where S = K - L , can also be partitioned into K smaller (LxL) 

sub-matrices, as follows: 



W = 



Eq(9) 



Each sub-matrix W;^ , for = {1 ... K} , is formed based on a different set of L rows of the matrix 
W . The concatenation of the K sub-matrices ,fork = {1 ... K}, would make up the matrix 
W. 



[0046] FIG. 3B graphically shows the partitioning of the vector H into K sub- vectors Hj^ , for 

k={l ... K} , and the partitioning of the matrix W into K sub-matrices Wj^ , for = {1 ... K} . 

[0047] It can be shown that the K sub-matrices are related to each other by the following 

relationship: 

Wjt=W,S, ,for A: = {2 ... K}, Eq (10) 

where is an (LxL) diagonal matrix that may be given as: 
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= diag 

which may be rewritten as: 
2^ =diag 



a-l)L (t-l)2L 

-/2>r — - — -/2>r 

I e N ^ N 



a-l)(L~l)L^ 



Eq(lla) 



1 e ^ 



... e 



(k-lHL-l)\ 



Eq(llb) 



As shown in equation (10), the K sub-matrices Wj^ , for A: = {1 ... K}, are related to each other, 
and the matrices W^^ , for ik = {2 ... K}, may each be derived based on the "base" sub-matrix 

[0048] It is observed that the diagonal elements of each matrix Zj^ , for /: = {1 ... K}, constitute a 

"generalized" column of a (TxT) DFT matrix whose elements are defined as shown in equation 
(2), except that N is replaced by T. For a (TxT) DFT matrix, the row index n and the column 
index m each run from 1 to T. However, for a generalized column of the (TxT) DFT matrix, 
the row index n can take on any integer value, and the elements of the generalized column would 
simply be repeated if and when the row index n exceeds T. In equation (1 lb), L may or may not 
be equal to T. The row index n for an (LxT) DFT matrix may then extend past the row 
dimension of the (TxT) DFT matrix if L>T, which would then result in the generalized 
column. 

[0049] A low-complexity least square (LCUS) method may be used to derive the least square 

channel impulse response estimate h'^ . For the LCLS method, equation (5) is first rewritten 
using the sub-vectors and the sub-matrices W^^ , for ^ = {1 ... K}^ as follows: 



Eq(12) 



Equation (12) may be expressed as the matrix product of an (LxL) matrix A and an (Lxl) 
vector B . The matrix A may be expressed as: 



A=|^|;wfW,j =(W''W)-' . 



Eq(13) 



Since the matrix A does not depend on the vector H, it can be computed off-line (i.e., pre- 
computed) and stored in a memory unit. 
[OOSO] The vector B may be expressed as: 
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Eq(14) 



Using the relationship for the sub-matrices Wj^ , for k = {l ... K}, shown in equation (10), the 
vector B may be rewritten as: 



1 0 

0 a', 



0 



... 0 
0 

0 



Eq(15) 



where w„ is the m-th column of the sub-matrix W, ; 

(il-lXm-l) 



a- =e 



-j2)T^ 



[0051] 



^ , for m = {1 ... L} ; and 

" denotes the conjugate transpose. 
Equation (15) may be simplified as follows: 



B = 



w. 



f ^ \ 



Eq(16) 



, As shown in equation (16), the vector B includes L inner products for the L entries of this 



vector. Each inner product is computed between the vector w" and the quantity 



obtain 



For each inner product, the quantity 



can be computed using 



one (TxT) DFT, as described below. 
[0052] A (TxT) DFT can be computed using a radix-2 fast Fourier transform (FFT), which 

requires Cy j^^y^j = (T • log2T)/2 complex multiplications. If the radix-2 FFT is used to compute 
for the vector B based on equation (16), then the number of required complex multiplications is 
C3 =L'[(T-log2T)/2 + L], where the second L (inside the bracket on the right side of the 

equation) is for the L complex multiplications needed for the inner product between w" and 
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^cx^M^k » the first L (on the right side of the equation but outside the bracket) is for the L 

inner products for the vector B . may also be expressed as Cg = + L - T • log2T/2 . The 
number of complex multiplications needed for the matrix multiply of the matrix A with the 
vector B is = L • L = . The total number of complex multiplications needed to compute 

the estimate h using the LCLS method and radix-2 FFT and based on equation (16) may then 
be expressed as: 

C«u.u,<ux2=CAB+CA=2L^+L.T.log,T/2 . Eq (17) 

[0053] For the example OFDM system described above, L = 16, T = 16, and K = 14 for the 

downlink. The total number of complex multiplications needed to compute the estimate h'^ 
using the brute-force method based on equation (5) is C^f = 16 • 16 • 14 = 3,584 . The total number 

of complex multiplications needed to compute the estimate h using the LCLS method with 
radix-2 FFT and based on equation (16) is Cj^^^^^ = 2-16^ +16-16-log2l6/2 = 1024. This 
represents a reduction of 71.42% in the number of complex multiplications required to compute 
the estimate h . 

[0054] The (TxT) DFT can also be computed using a radix-4 FFT, which requires 

Ct^^x4 =((T/4-l)/(T/2)(T-log2T) complex multiplications. The total number of complex 

multiplications needed to compute the estimate h using the LCl^ method and radix-4 FFT and 
based on equation (16) is C^^^_^^^^ = 16 • 16 + 16 • [(3/8)(16 • logjie) + 16] = 896 . This represents 
a reduction of 75% in the number of complex multiplications required to compute the estimate 
h . Table-1 lists-the number of complex multiplications required to compute the channel impulse 

response estimate h using (1) the brute-force method and (2) the LCLS method with radix-2 

and radix-4 FFFs. also shows the percentage savings achieved by the LCLS method over the 
brute-force method. 
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Table 1 





Number of complex multiplications 


Savings 
(%) 


Brute-force (Cbf) 


LCLS (Ctotal) 


Radix-2 FFT 


3,584 


1,024 


71.42 


Radix-4 FFT 


3.584 


896 


75.00 



[0055] FIG. 4 is a flow diagram of an embodiment of a process 400 for deriving a least square 

channel impulse response estimate using the low complexity least square method described 
above. In the following description, the inputs to the process are as follows: 

• Number of taps for the channel impulse response: L; I 

• Number of total subbands: N = L • T ; 

• Number of designated subbands: S = L • K ; and 

• Initial channel frequency response estimate H with channel gains for the S designated 



subbands: H = [Hp H 



F+l 



H 



F+LK- 



The output of the process is the least square channel impulse response estimate, 

h =[h, h, ... 

[0056] Initially, the S entries of the (Sxl) vector H are arranged into a (TxL) matrix Htxl 

(step 412), as follows: 

H 



HtxL "~ 



H 



F+i 



H 



F+L 



H 



F+L+1 



H 



F+(L-1)K 

0 



F+(L-1)K+1 

0 
0 



H 



H 



^ F+L-l 
F+2L-1 

F+LK-1 

0 
0 



Eq(18) 



- 0 

As shown in equation (18), the S entries of the vector H are written row-wise into the matrix 
H^j^L, starting in the first row and going from left to right. Each row of the matrix H^xl 
includes L consecutive entries of the vector H . The matrix Hj^l thus effectively partitions the 
vector H into K (Lxl) sub- vectors Hj^ , for /: = {!... K}, where each sub- vector H^^ 
corresponds to one row of the matrix H^xl • The matrix H^xl includes N entries for the N total 
subbands. Since S is typically less than N, only the first K rows of the matrix H^xl include non- 
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zero values ifrom the vector H and the last entries in the matrix H^xl ^ filled with 

zeros, as shown in equation (18). 
[0057] An (LxT) DFT matrix W^^t is next formed (step 414). The (n,m) -th entry, wj^, of 

the matrix Wlxt is defined as: 

.(n-lXw-l) 



wj^ = e T , for n = {1 ... L} and m = {1 ... T} . Eq (19) 



Each column of the matrix Wl^t corresponds to a generalized column of a (TxT) DFT matrix. 
The m-th column of the matrix Wlxt ^^us includes L entries that correspond to the diagonal 
elements of the matrix 1,^ shown in equation (lib), where k=m for m = 1 through K. Since the 
index k for the matrices from 1 through K but the index m for the columns of the matrix 

Wlxt ru"S from 1 through T, and since K < T , not all columns of the matrix Wlxt ^^^y 

used, as described below. 

^^^^ A 
[0058] FIG. 3C graphically shows the relationship between the matrices Wlxt Htxl • Th® 

first K rows of the matrix Htxl correspond to the K sub-vectors H^^ 'for k = {1 ... K} . The T 

columns of the matrix Wt^-r are generalized columns of a (TxT) DFT matrix. Each row of the 

matrix Wlxt is a "normal" row of the (T xT)^ DFT matrix. 

[0059] T-point DFTs of the columns of the matrix H^xl ^® then computed using the matrix ^ 

Wlxt (step 416), The DFTs may be expressed as: 



Glxl — Wlxt Htxl 



i. 

i2 



Eq(20) 



where. g._,.for..m..=.{l .L] , is.an_(Lx l.)-rowr_Yector-for-the-m---th-row-of-the matrix-GLx,;^- 



[0060] Each row-vector g includes L entries, where each entry is obtained based on a T-point 

— m 

\ 

DFT of one row of the matrix Wlxt one column of the matrix HtxL' ^s shown in FIG. 3C. 
Equation (20) essentially performs the computation for the L sunmiations shown in equation 
(16); such that 

i.=Z<H* ,form = {l,.. L}, Eq(21) 
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[0061] The matrix Wi^t includes T columns for T generalized rows of a (TxT) DFT matrix. 

However, only the first K columns of the matrix W^xt are used for the K matrices Zj^ , for 

k = [l ... K}. The last (T-K) columns of W^^t. are not used, since these columns are 

multiplied with the last (T — K) rows of zeros in the matrix Hj^l • 
[0062] Each of the L entries of the vector B is then obtained by computing an inner product 

between conjugate transpose of a vector w„ and a corresponding row- vector (step 418). 

This inner product may be expressed as: 

ft.=5i::gl=g„w: ,form = {l...L}, Eq (22) 

where is the m-th column of the (LxL) sub-matrix W, , and " * " denotes a conjugate. The 

sub-matrix W, is defined such that the (n,m) -th entry, , is given as: 

wl^=e N for ^ = L} and m = {l ... L}. Eq (23) 

The result of step 418 is the vector B = [b^ ^lI^ • 

[0063] The (LxL) matrix A may be pre-computed as shown in equation (13) and stored in a 

memory unit (step 420). The least square channel impulse response estimate h may then be 

computed by performing a matrix multiply of the matrix A with the vector B (step 422). This 
matrix multiplication may be expressed as: 

h'' = AB . Eq(24) 

[0064] FIG. 5 is a flow diagram of another embodiment of a process 500 for deriving a least 

square channel impulse response estimate using the low complexity least square method. 

[0065] Initially, an intermediate vector is derived based on (1) K sub-vectors of a vector for a 

first channel estimate and (2) at least two DFT sub-matrices for a DFT matrix (step 512). The 

intermediate vector may be B , the K sub- vectors may be Hj^ , for k={l ... K}, the vector for 

the first channel estimate may be H, the first channel estimate may be the initial channel 

frequency response estimate, the at least two DFT sub-matrices may be W^^ , for /: = {1 ... K}, 

and the DFT matrix may be W. The intermediate vector B may then be obtained by (1) 
performing a matrix multiply of each of the K sub-vectors with a corresponding one of the K 
DFT sub-matrices to obtain a corresponding intermediate sub- vector Wj^Hj^, and (2) 
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accumulating K intermediate sub- vectors Wf H^, for A: = {1 ... K}, to obtain the intermediate 
vector B , as shown in equation (14). 
[0066] Alternatively, the at least two DFT sub-matrices may be Wi^^. and Wj. The 

intermediate vector B may then be obtained by (1) computing DFTs of a first matrix H^xl, 
formed based on the vector H for the first channel estimate, to provide a second matrix G^xl » 
and (2) computing inner products between the colunms of the base DFT sub-matrix Wj and the 
rows of the second matrix G^xl obtain the intermediate vector B , as shown in equations (21) 
and (22). 

[0067] An intermediate matrix is then obtained, which is derived for the DFT matrix 

corresponding to the vector for the initial frequency response estimate (step 514). The 
intermediate matrix may be the matrix A, which may be derived as shown in equation (13). 
Again, the matrix A may be pre-computed, stored in a memory unit, and retrieved when 
needed. 

[0068] A second response estimate is then derived based on the intermediate vector and the 

intermediate matrix (step 516). The second response estimate may be a least square channel 
impulse response estimate. 

[0069] In the above description, the structure of the matrix W is exploited to greatly reduce the 

complexity of the derivation of the least square channel impulse response estimate. The reduced 
complexity channel estimation techniques described herein may also be used to derive other 
channel estimates. For example, these techniques may possibly be used to derive the enhanced 
channel frequency response estimate shown in equation (6). In general, these techniques may be 
used for any problem where multiplication by a sub-matrix of the DFT matrix is involved. 
However, the gains achieved by these techniques may be dependent on the set-up of the 

problemr ~ 

[0070] As noted above, the channel estimation techniques described herein may be used for any 

communication system with multiple subbands, such as OFDM systems. Moreover, these 
techniques may be used for multiple-input multiple-output (MIMO) systems that employ 
multiple (N^. ) transmit antennas and multiple ( Nr ) receive antennas for data transmission. For 
a MIMO system that utilizes OFDM, the response of a wireless MIMO channel may be given as 
H(/:), for k={l ... N}. Each matrix H(k) is an (N^xN^) matrix with entries Z/.^, for 

I = {1 ... Nr } and j = {1 ... N^.} , where //. ^. is the channel gain between the j-transmit antenna 
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and the i-Xh receive antenna. The techniques described herein may be used to derive the channel 
response of each transmit/receive antenna pair. 

[0071] FIG. 6 is a block diagram of an embodiment of an access point 600 and a terminal 650, 

which are capable of deriving the channel estimate using the techniques described herein. 

[0072] On the downlink, at access point 600, traffic data is provided to a TX data processor 610, 

which formats, codes, and interleaves the traffic data to provide coded data. An OFDM 
modulator 620 then receives and processes the coded data and pilot symbols to provide a stream 
of OFDM symbols. The processing by OFDM modulator 620 may include (1) symbol mapping 
the coded data to form modulation symbols, (2) multiplexing the modulation symbols with pilot 
symbols, (3) transforming the modulation symbols and pilot symbols to obtain transformed 
symbols, and (4) appending a cyclic prefix to each transformed symbol to form a corresponding 
OFDM symbol. For the downlink, the pilot symbols may be multiplexed with the modulation 
symbols using, for example, time division multiplexing (TDM). For TDM, the pilot symbols 
and modulation symbols are transmitted on different time slots. The pilot symbols may be 
transmitted on designated subbands, where may include all or a subset of the M usable 
subbands. 

[0073] A transmitter unit (TMTR) 622 then receives and converts the stream of OFDM symbols 

into one or more ^analog signals and further conditions (e.g., amplifies, filters, and frequency 
upconverts) the analog signals to generate a downlink modulated signal suitable for transmission 
over the wireless channel. The downlink modulated signal is then transmitted via an antenna 
624 to the terminals. 

[0074] At terminal 650, the downlink modulated signal is received by antenna 652 and provided 

to a receiver unit (RCVR) 654. Receiver unit 654 conditions (e.g., filters, amplifies, and 
frequency downconverts) the received signal and digitizes the conditioned signal to provide 
samples. An OFDM demodulator 656 then removes the cyclic prefix appended to each OFDM 
symbol, transforms each recovered transformed symbol using an FFT, and demodulates the 
recovered modulation symbols to provide demodulated data. An RX data processor 658 then 
decodes the demodulated data to recover the transnwtted traffic data. The processing by OFDM 
demodulator 656 and RX data processor 658 is complementary to that performed by OFDM 
modulator 620 and TX data processor 610, respectively, at access point 600. 

[0075] OFDM demodulator 656 may further determine an initial frequency response estimate 

for the downlink channel, or provide the received pilot symbols that may be used to derive 
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. A processor 670 receives (or equivalent information) and may derive a least square 

impulse response estimate of the wireless channel based on H^j,, and using the low 

complexity least square method described above. Processor 670 may further obtain an enhanced 

^ is ^is 

frequency response estimate for the downlink channel based on h^. The enhanced 

^ Is 

estimate H.^ may thereafter be used for uplink data transmission and/or sent back to the access 
point for use for downlink data transmission. 
[0076] On the uplink, traffic data is processed by a TX data processor 682 and provided to an 

OFDM modulator 684, which also receives pilot symbols. OFDM modulator 684 may then 
process the coded data and pilot symbols similar to that described for OI^M modulator 620. 
For the uplink, the pilot symbols may also be multiplexed with the modulation symbols using 
TDM. Moreover, the pilot symbols may be transmitted on only S^p . subbands that have been 

assigned to terminal 650 for pilot transmission. 

[0077] A transmitter unit 686 then receives and processes the stream of OFDM, symbols to 

generate an uplink modulated signal suitable for transmission over the wireless channel. The 
modulated signal is then transmitted via an antenna 632 to the access point. 

[0078] At access point 600, the uplink modulated signal is processed by a receiver unit 642 to 

provide samples. These samples are then processed by an OFDM demodulator 644 to provide 
demodulated data, which are further processed by an RX data processor 646 to recover the 
transmitted traffic data. OFDM demodulator 644 may determine the initial frequency response 

estimate H^p , for the uplink channel for each active terminal or provide the received pilot 

synibols that may be used to obtain H„p , . A processor 630 receives H^p , (or equivalent 
information) for each active terminal, determines the least square channel impulse response 
estimate h^p , for the active terminal based on Hup , and using the low complexity least square 

^ Is 

method, and further obtains the enhanced channel frequency response estimate Hyp based on 

h„p ; . The enhanced estimate H^j may thereafter be used for downlink data transmission to the 

terminal and/or sent back to the terminal for use for uplink data transmission. 
[0079] Processors 630 and 670 direct the operation at the access point and terminal, respectively. 

Memory units 632 and 672 provide storage for program codes and data used by controllers 630 
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and 670, respectively. Processors 630 and 670 may be designed to perform the computation 
described above to derive estimates of the uplink and downlink channels, respectively. 

[0080] The reduced complexity channel estimation techniques described herein may be 

implemented by various means. For example, these techniques may be implemented in 
hardware, software, or a combination thereof. For a hardware implementation, the elements 
used to implement any one or a combination of the techniques may be implemented within one 
or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital 
signal processing devices (DSPDs), programmable logic devices (PLJDs), field progranmiable 
gate arrays (FPGAs), processors, controllers, micro-controllers, microprocessors, other 
electronic units designed to perform the functions described herein, or a combination thereof. 

[0081] For a software implementation, the channel estimation techniques may be implemented 

with modules (e.g., procedures, functions, and so on) that perform the functions described 
herein. The software codes may be stored in a memory unit (e.g., memory units 632 or 672 in 
FIG. 6) and executed by a processor (e.g., processor 630 or 670). The memory unit may be 
implemented within the processor or external to the processor, in which case it can be 
conununicatively coupled to the processor via various means as is known in the art. 

[0082] The previous description of the disclosed embodiments is provided to enable any person 

skilled in the art to make or use the present invention. Various modifications to these 
embodiments will be readily apparent to those skilled in the art, and' the generic principles 
defined herein may be applied to other embodiments without departing from the spirit or scope 
of the invention. Thus, the present invention is not intended to be limited to the embodiments 
shown herein but is to be accorded the widest scope consistent with the principles and novel 
features disclosed herein. 



